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Using cluster perturbation theory, we explain the origin of the strongly dispersive feature found 
at high binding energy in the spectral function of the Hubbard model. By comparing the Hubbard 
and t—J—3s model spectra, we show that this dispersion does not originate from either coupling to 
spin fluctuations (oc J) or the free hopping (oct). Instead, it should be attributed to a long-range, 
correlated hopping oc /U, which allows an effectively free motion of the hole within the same 
antiferromagnetic sublattice. This origin explains both the formation of the high energy anomaly 
in the single-particle spectrum and the sensitivity of the high binding energy dispersion to the 
next-nearest-neighbor hopping t'. 

PACS numbers: 71.10.Fd, 74.72.Cj, 79.60.-1 


INTRODUCTION 

High-temperature superconductivity in the cuprate 
oxides has attracted significant attention over the past 
thirty years. However, the precise origin of this phe¬ 
nomenon is not well understood due to the complex 
physics even in a minimal model used to describe the 
correlated nature of the electrons [1-3] : the 2D Hubbard 
model. It is often assumed that a first step in under¬ 
standing the physics of this model is to study its spectral 
properties [4-17] in the simple undoped limit. 

The spectral function of the undoped 2D Hubbard 
model when the free electron bandwidth W is comparable 
to the Hubbard interaction U consists of two prominent 
features in the band structure, cf. Fig. 1. At low binding 
energies (LBE), a sharply defined quasiparticle-like exci¬ 
tation disperses downward from (7r/2,7r/2) reaching an 
energy of the order of spin exchange J = At'^/U near the 
F-point (0,0). This quasiparticle, often labeled the spin 
polaron (SP) [18], represents a hole heavily dressed by 
spin excitations from the antiferromagnetic (AF) ground 
state with its bandwidth no longer governed by the free 
electron hopping t but by the spin exchange J [9, 18-23]. 

At higher binding energies (HBE), another feature 
becomes prominent approaching the F-point. The delin¬ 
eation of the LBE spin polaron from the HBE feature has 
been widely associated with the “high energy anomaly” 
or “waterfall” seen in a number of cuprate photoemission 
results [11, 21-40]. While the spin polaron is well 
understood, the characterization remains poor for the 
feature at higher energies. One may naively expect that 
a broad, non-dispersive Hubbard band should appear at 



FIG. 1: (color online) Left panel: spectral function A(k,a;) 
of the undoped Hubbard model with U = 8t and t' = 0 
calculated using CPT on a 4 x 4 cluster for the high symmetry 
directions of the 2D Brillouin zone. The blue dotted (green 
dashed) box guides the eye for the LBE (HBE) features. Right 
panel: the high symmetry directions of the 2D Brillouin zone. 


the F-point. However, as shown in Fig. 1, it is clear 
that there is a sharp dispersion. Previous studies have 
attributed this feature to scenarios such as spin-charge 
separation [16, 17, 34, 35, 41-45] or a weak-coupling spin- 
density wave [9, 10, 46]. However, these interpretations 
remain controversial. 

The aim of this paper is to understand the nature of the 
HBE feature and what separates it from the SP. Using 
cluster perturbation theory we examine both the Hub¬ 
bard, t—J and t—J—3s models to provide an in-depth 
examination of what controls the quasiparticle dispersion 
and intensity at low and high energies. We demonstrate 
that the HBE dispersion is mainly set by a correlated 
effective hopping /U rather than the free hopping t. 
While the SP dispersion is set by spin exchange J due to 
the strong coupling of the mobile hole to spin excitations. 















2 



FIG. 2: (color online) Spectral function ^(k,a;) for the 
undoped Hubbard model with U = 8t and (a) t' = —0.3t 
and (b) t' = 0.3t calculated using CPT on a 4x4 cluster. 


the HBE dispersion is determined by hopping within the 
same AF sublattice through long-range (so-called 3-site) 
correlated hopping. Therefore, one naturally expects a 
distinct transition between them, providing a physical 
picture for the high energy anomaly. Moreover, as the 
uncorrelated next-nearest-neighbor (nnn) hopping t' also 
allows for hopping within the same AF sublattice, we 
obtain a natural explanation for the sensitivity of the 
Hubbard model spectrum at high energies to changes in 
t'. 



FIG. 3: (color online) (a) Upper panel: schematic cartoon 
showing the 3-site correlated hopping elements oc t^/t/ in an 
exchange process via virtual double occupancy oc U\ lower 
panel: various 3-site paths in the AF background. Arrows 
(circles) denote electron spins (holes), (b-c) Spectral function 
A(k,a;) for the undoped (b) t—J and (c) t—J—3s model with 
J = 0.5t calculated using GPT on a 4x4 cluster, (d) The 
same as (c) but with the spectral weight normalized in the 
Hubbard-like way, see text. 


MODELS AND METHODS 

The 2D Hubbard model Hamiltonian Hh is 
'M-h ^ ^ “b h.C^ 

( 1 ) 

where (cia) and Uicr denotes the creation (annihila¬ 
tion) and density operator at site i of spin a; ty is the 
hopping, restricted here to nearest-neighbors (nn) =t 
and nnn t((y)) =t'. The Hubbard spectral function has 
been calculated using various numerical methods, e.g. 
exact diagonalization (ED) [5-7], quantum Monte Carlo 
(QMC) [8-11] or cluster perturbation theory (CPT) 
based on ED [12-17]. Although CPT is an approximate 
method, we believe that for the purpose of this paper it is 
most suitable, as it can produce continuous momentum 
resolution at zero-temperature. 

It is well known that t' qualitatively changes the 
quasiparticle dispersion and mimics the differences be¬ 
tween hole- and electron-doping in cuprates found using 


Hamiltonians which incorporate higher energy degrees of 
freedom [11, 21-35], cf. Fig. 2. A positive (negative) t' 
is found to increase (decrease) the high-energy anomaly 
which shows that we can indeed identify two distinct 
features (LBE and HBE) in the spectrum. Moreover, 
this comparison also shows that the HBE dispersion is 
far more affected by t' than is the SP and that changing 
t' does not introduce new qualitative features in the 
spectral function. This suggests that the dispersion 
should somehow be related to an effective nnn hopping 
which must be present even when t' = 0. 

While this effective long-range hopping is not explicitly 
present in the Hubbard model when T = 0, it does appear 
in a t/U expansion to lowest order [47-50] as the so- 
called 3-site term, see Fig. 3(a) and Refs. 49-64. This 
defines the t—J—3s model with the Hamiltonian given 
by Ht-j-ss = JLt-j + 'Hzs '■ 

Ut-j = -^ E (cj,,G<T + /i.c.)+J 

(ij> 

T~1-3s = ^ E! , (2) 































3 


where S; • Sj = 5'fS'? + i (^S^Sr + with S? = 

(n;^ — 714)72 and = (Sr)^ = Cj^c 4 . The constrained 
fermionic operators acting in the Hilbert space without 
double occupancies are defined as C;^ = ni^)- Note 

that only when = J does this t—J—3s model follow 
from a perturbative expansion of the Hubbard model, 
Eq. (1). 


SPECTRAL FEATURES 

To verify that the 3-site terms are responsible for the 
onset of the HBE dispersion, we calculate the spectral 
function of both the t—J and t—J—Zs model using CPT, 
see Figs. 3(b) and (c) respectively [65]. At low energies 
the t—J and t—J—Zs model spectra are quite similar, 
see Refs. 53, 66 , and qualitatively reproduce the SP 
dispersion in the Hubbard spectrum. The qualitative 
agreement extends to the HBE only in the t—J—Zs 
model, showing explicitly that the 3-site terms [52- 
57, 61-64] indeed play a crucial role in the development 
of HBE dispersion. We note that the spectral weight 
of the t—J—Zs and Hubbard models remains quite 
different which can be attributed to different sum rules. 
Artificially introducing equivalent spectral weight sum 
rules between the two models [see Fig. 3(d)] produces 
an even more qualitative and quantitative agreement 
between the two spectra (see Appendix for a detailed 
discussion). 

With a correlation established between the HBE dis¬ 
persion and the 3-site terms in a t/U expansion of the 
Hubbard model, a natural question is whether or not the 
dispersion is real and, if so, how does one understand the 
underlying mechanism which dominates this dispersion. 
We postulate that the 3-site terms open an additional 
channel through which a hole can “freely” propagate 
within the same AF sublattice. This motion does not 
disturb the underlying AF order, as shown intuitively 
in Fig. 4(a), and we call this mechanism an intra- 
AF-sublattice effective hopping (lASH). The dispersion 
relation is given by 

^QQs2kx+cos2ky+4:COskx cosfc^], (3) 

which can be deduced from the Hsg term assuming a 
perfect AF background. 

To test this postulate and provide an answer to these 
questions, we calculate A(k, w) for the t—J—Zs model 
using a more analytical technique - the self-consistent 
Born approximation (SCBA) [18, 21, 53, 67]. This 
approximate technique is not limited to small cluster 
sizes [18] and can be viewed as complementary to CPT. 
To emphasize the role of the 3-site terms, we calculate 
the t—J—Zs model using SCBA for various values of 
Jss as shown in Figs. 4(b-f), where the physical value 
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M 


FIG. 4: (color online) (a) Idealized propagation via lASH 
(dashed lines) triggered by e.g. the 3-site terms: the hole 
motion does not disturb the AF background, (b-f) Spectral 
functions A(k, ca) for various 3-site term values in t—J—Zs 
models using SCBA on a 40x40 cluster (J = 0.5t and 
broadening <5 = 0.15t), with Jss = 0, J/2, J, 3J/2 and 4J, 
respectively. [The canonical parameter Jss = J is highlighted 
by the red box, panel (d)]. Dotted lines indicate fitted 
SP dispersion, while dashed lines indicate lASH dispersion 
Eq. (3). 


of J 3 S obtained from the perturbative expansion of 
Hubbard model is Jss = J [49-60] [see Fig. 4(d)]. 
First, with increasing J^g a continuous spectral weight 
transfer occurs from the top of the SP to higher energies, 
proceeding from the P and M points [21-23, 53]. Next, 
comparing Fig. 4(d) with Fig. 3(c), we conclude that the 
HBE dispersion in both CPT and SCBA shows a very 
specific shape around the T point with the “proper” value 
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of the 3-site terms. This suggests that such a dispersion 
is “real” and not an artifact of the CPT method which 
artificially may have enhanced the dispersion due to the 
approximation. Neglecting the 3-site terms [Fig. 4(a)] 
suppresses the signatures of dispersion at high energy [68] 
with the only prominent dispersion associated with the 
SP (see Appendix and Refs. 18, 21, 53 for a detailed 
discussion). As a single dispersive branch is visible for 
an artificially large value of J^s [see Fig. 4(f)], there is 
likely an onset of the dispersion due to the 3-site terms. 

The HBE dispersion always is well-approximated by 
Eq. (3) [cf. Fig. 4(b-f)]: for smaller J^s values this is 
visible around the F and M points and the region of 
agreement grows across the Brillouin zone for larger 
Note that this dispersion relation is deduced from Tiss 
assuming a perfect AF background and neglecting the 
coupling to quantum spin fluctuations. This suggests 
that in the corresponding parts of the Brillouin zone 
the hole indeed can be mobile, unencumbered by spin 
fluctuations, as illustrated by the cartoon in Fig. 4(a). 
This is qualitatively different from the situation in the 
t—J^—2>s model, as discussed in the Appendix and 
Refs. 62, 63. 


CONCLUSIONS AND DISCUSSION 

In summary, we have presented a detailed understand¬ 
ing of the spectral function of the undoped Hubbard 
model, identifying distinct origins for the low and high 
binding energy dispersions. Comparing the spectral 
functions of the Hubbard, t—J and t—J—3s models, 
we have argued that the dispersion at high binding en¬ 
ergy originates from an effective intra-antiferromagnetic- 
sublattice hopping, primarily due to the so-called 3- 
site terms. We found that this dispersion relation is 
not renormalized by coupling to spin fluctuations which 
suggests that holes at these energy scales can propagate 
more-or-less freely on the same sublattice. This distinc¬ 
tion between the physical mechanisms which give rise to 
the dispersion at both low and high energies, and thus 
the character of the associate quasiparticles, unveils a 
rather profound origin for the high energy anomaly. The 
result also has important consequences for the Hubbard 
model with uncorrelated nnn hopping t', explaining its 
crucial impact on the dispersion primarily at high binding 
energy. 

Given the clear dichotomy between the features of 
the spectral function at low and high energy and the 
correlated, but effectively free, hopping associated with 
the 3-site terms, we might expect distinct effects on 
emergent behavior, e.g. antiferromagnetism, charge or¬ 
der, superconductivity, or some combination as found in 
the cuprates. In particular, since t' or nnn hopping has 
been linked to Tc [69, 70], this may further indicate a 
causal connection between the HBE feature and material 


dependence of Tc- Moreover, the fact that a doped 
hole can effectively decouple from the antiferromagnetic 
background and move without disrupting the spin order 
in part of the Brillouin zone is consistent with the re¬ 
cently observed persistence of magnetic excitations upon 
doping [71-76]. Since none of the prominent dispersing 
features in the half-filled Hubbard model are due to 
the free hopping t, a natural follow-up question would 
be when and how one might expect the uncorrelated 
hopping t to become important upon doping? Since there 
should be a crossover between the Mott insulator near 
half-filling and a weakly-correlated metal at substantial 
doping, the answer to this question could shed some light 
on our understanding of superconductivity which appears 
near optimal doping in the cuprates. 
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Details of the CPT Calculations 


Evaluation of the spectral function The single particle 
spectral function of any quantum mechanical Hamilto¬ 
nian H is defined as: 


A(k, ui) = -Im 

TT 




(Gjc; 


1 


k(7 


uj + T-L — Ec -f iS 


TyCkrrlG), 


( 4 ) 

where (cko-) denotes the electron creation (annihi¬ 
lation) operator with momentum k and spin cr, jG) is 
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the ground state with energy Eg, and 5 is a Lorentzian 
broadening. 

On a finite cluster, the Hamiltonian can be split 
into H = He + Hint, where He contains the (open¬ 
boundary) intra-cluster terms while Hint contains the 
inter-cluster hopping terms. The CPT method evaluates 
He exactly using exact diagonalization, and treats Hint 
perturbatively. The corrected Green function can be 
expressed as 


g = 


Ge 




(5) 


where Gc is the cluster Green function (evaluated from 
He)- In the long-wavelength limit 


0cpr(k,a;) = ^ (g) 

a,b 


where a, b are intra-cluster site indices. Under the CPT 
method, the spectral function is evaluated by 

A{k,uj) = -C/c'PT(k,w). (7) 

TT 

Justification of using CPT Although CPT is not 
strictly applicable to t—J-like models because the two- 
particle vertex (spin exchange on the boundary) is not 
well defined, the t vertex should be dominant in an intra¬ 
cluster t/U expansion of the Hubbard model. Further¬ 
more, for these quantum cluster approaches the key to 
the quality of the results lies in the cluster sizes, instead 
of higher-order corrections to the CPT formula. We 
also have studied the t—J—3s model on various cluster 
sizes such as 20 and 24 sites (not shown) and find that 
the dispersion saturates with increasing cluster size. In 
addition, one can assume that the agreement between 
the Hubbard and t—J—3s model spectra (Figs. 1 and 3) 
cannot be a coincidence (compared over different cluster 
sizes), which provides an a posteriori justification that 
the CPT method with only first order t-vertex can be 
acceptable for the t—J models. 

“Normalization” of the t—J—3s spectra Figure Al(a) 
shows the momentum distribution of integrated spectral 
weight (electron density) nt in the lower Hubbard band 
(below Ep). This is different in the t—J—3s model, 
where due to the projection of the upper Hubbard 
band (forbidden double-occupancy), the spectrum has a 
uniform nk distribution. Thus, to partially correct for 
this difference and enhance comparisons, the Fig. 3(d) is 
obtained by dividing the spectrum found in Fig. 3(c) by 
the Uk distribution for the Hubbard model [Fig. Al(a)] 
at each momentum. 

Quantitative comparison of CPT results In what fol¬ 
lows we show in a quantitative way that: (i) the CPT 
results for the t—J—3s model are closer to those of the 
Hubbard model than the CPT results for the “bare” t—J 
model, and (ii) the CPT results for the t—J—3s model, 




FIG. Al: (color online) (a) The electron density rik used 
for normalization in Fig. 3. (b) The momentum-integrated 
correlation function ^{lj) [Eq. (8)] calcnlated for the CPT 
A(k,aj) resnlts for: the t—J and Hubbard models (dot- 
dashed lines), t—J—3s and Hubbard models (dashed lines) 
and t—J—3s models normalized in the Hubbard-like way [or 
{t—J—3s)/N] and Hubbard models (solid lines). 


normalized in a Hubbard-like way, qualitatively repro¬ 
duce the Hubbard model spectral function extremely well 
{i.e. better than unrenormalized t—J—3s results). 

We define the (momentum-integrated) correlation 
function of two spectral functions A(“^(k, w) and 
A(«(k,w): 




/A(“)(k,a;)A(^)(k,a;)d2k 
|A(“)(k,w)|2d2k/ |A(/3)(k,w)|2d2k 


• ( 8 ) 


Fig. Al(b) shows the correlation function between 
the Hubbard spectrum, with t—J , t—J—3s , and 
normalized t—J—3s spectrum calculated via CPT. One 
concludes that: (i) indeed the “normalized” t—J—3s 
model reproduces extremely well the spectrum of the 
Hubbard model also on a semi-quantitative level, (ii) 
the “unnormalized” t—J—3s model provides a marginally 
“worse” comparison than the “normalized” t—J—3s 
model, and (hi) the “bare” t—J model certainly is unable 
to capture the features, especially at high binding energy 
(~ oj — p < —4t), of the Hubbard model. 


Details of the SCBA Calculations 

To calculate the spectral function A(k, w) for the 
undoped t—J—3s model, we first assume that the ground 
state has (i) a broken spin rotational symmetry and (ii) 
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long range AF order. These two assumptions allow us to 
map the t—J—3s model onto the so-called spin-polaron 
Hamiltonian, cf. Eq. (2.15) in Ref. 53. Next the spectral 
function for the spin polaron Hamiltonian is calculated by 
using the SCBA: we calculate the self-energy by summing 
all noncrossing diagrams to the infinite order and self- 
consistently evaluate the single particle Green’s function 
on 40x40 lattice sites. 


t-J^ 



FIG. Cl: (color online) Comparison of A(k, a;) for t—Jz—3s 
model calculated using SCBA for various strength of the 3- 
site term (a) Jss = 0, (b) Js* = J/2, (c) Jss = J, (d) 
Jss = 2J and (e) Jsa = 4J, see text for further details. The 
dashed line indicates the lASH dispersion following from the 
inclusion of the 3-site terms in the model. 


We also present here the functional form of the SP dis¬ 


persion relations following the considerations of Ref. 18: 

= xi + a; 2 (cos + cos ky)"^ 

-l-xa [cos(2fca;) + cos{2ky)], (Bl) 

where xi, X 2 , and 0:3 are parameters which depend 
on the ratio of J/t in the t—J-like model and which 
follow from fitting the above dispersion relation to the 
dispersion of the upper edge of the continuum in the 
SCBA calculated spectra. In the case of the t—J model 
we obtain xi = —1.5t, X 2 = 0.18t, and 0:3 = 0.035t (note 
a slightly different constant shift xi with respect to the 
one adopted in Ref. 18 due to a different system size in the 
SCBA calculations). In the case of SCBA calculations of 
the t—J—3s model we obtain xi = —1.6t, X 2 = 0.33t, and 
0:3 = 0.035t. Note that the difference between the fitting 
parameters for t—J and t—J—3s models suggests that 
the SP dispersion relation is influenced by the presence 
of the 3-site terms. This suggests that the hole forming 
the SP moves not only as a result of coupling to the spin 
fluctuations but also via the 3-site terms. 


Influence of Quantum Fluctuation on the lASH 

We calculate the spectral function A(k, w) using SCBA 
for the t—Jz—3s model with varying strength of the 3-site 
term where 


(ij> 

(iJ).(iJO 

Compared with Fig. 4, the result presented in Fig. Cl 
shows that when quantum spin fluctuations are switched 
off the lASH shows up only for unrealistically large 
strength of the 3-site hopping term. As shown in Fig. Cl, 
the spectrum of the t—Jz—3s model shows almost no sign 
of an unrenormalized free dispersion which should result 
from the inclusion of the 3-site terms. Instead, each of 
the polaronic-like states, visible in the ladder spectrum, 
acquires a renormalized dispersion due to the onset of the 
3-site terms. In this case the lASH and the polaronic- 
like propagation are not independent of one another, as 
already explained in detail some time ago in Refs. 62, 63. 
It is then only in the case of an unrealistically large value 
of the 3-site terms that the unrenormalized dispersion 
appears in the spectral function of the t—Jz—3s model. 

This discrepancy between the t—Jz—3s and t—J—3s 
models is in contrast with the interesting result of Ref. 77. 
In that case the reported irrelevance of quantum spin 
fluctuations on the hole motion in a charge transfer 
insulator (as governed by the three-band model) concerns 
solely the effective SP dispersion and is due to the 
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rather similar dispersion relation of the quantum-spin- 
fluctuation-mediated hopping and the lASH. 
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